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Abstract. Spiking Neural P systems, SNP systems for short, are bio- 
logically inspired computing devices based on how neurons perform com- 
putations. SNP systems use only one type of symbol, the spike, in the 
computations. Information is encoded in the time differences of spikes or 
the multiplicity of spikes produced at certain times. SNP systems with 
delays (associated with rules) and those without delays are two of sev- 
eral Turing complete SNP system variants in literature. In this work we 
investigate how restricted forms of SNP systems with delays can be sim- 
ulated by SNP systems without delays. We show the simulations for the 
following spike routing constructs: sequential, iteration, join, and split. 
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1 Introduction 

Membrane computing abstracts computational ideas from the way biological 
cells perform information processing. The models used in Membrane computing 
are known as P systems. Spiking Neural P systems, or SNP systems for short, 
are a class of P systems that incorporate the neural-like P systems which only 
use indistinguishable electric signals or spikes. SNP systems were introduced in 
[5] and were proven in [5] and [3] to be Turing-complete. Because SNP systems 
only use one kind of spike symbol a, the output of the computation done by the 
system is based on the time interval between the first spike and the second spike 
of a designated output neuron. This model represents the fact that spikes in a 
biological neural system are almost identical from an electrical point of view. 
The time intervals between these signal spikes are crucial to the computations 
performed by neurons. 

As with the more common results in Membrane computing, SNP systems 
have been used to efficiently solve hard problems such as the SAT problem 

1 a relatively new field of Unconventional and Natural computing introduced in 1998 
by Gheorghe Paun, see e.g. [10] 



using exponential workspace [9]. Other numerous results also focus on using 
SNP systems as acceptors, generators, and transducers as in [7] and [T]. SNP 
systems were also used in order to perform arithmetic operations where the 
operands are encoded in their binary formats [5J. It has been proven in [3] that 
delays in applying rules in the neurons of an SNP system are not required for 
SNP systems to be Turing complete. However SNP systems with delays might 
still prove to be useful for modeling purposes. In [13 j SNP systems without 
delays were represented as matrices and the computations starting from an initial 
configuration can be performed using linear algebra operations. By simulating 
SNP systems with delays using SNP systems without delays, the work in (T3] 
can be used to further study SNP systemsj^] 

In this work we investigate how SNP systems with delays can be simulated 
by SNP systems without delays. In particular we focus on a special type of SNP 
systems that consist of a source neuron and a sink neuron. The initial config- 
uration consists of a spike in the source neuron only and a final configuration 
where only the sink neuron has a spike. The spike from the source neuron to 
the sink neuron is routed through other intermediary neurons using four routing 
primitives: sequential, iteration, joins, and splits. 

The contents of this paper are organized as follows: Section [2] provides def- 
initions and notations that will be used in this work. Section [3] presents our 
results. We end this paper in Section [4] with some final notes, conjectures and 
open problems. 



2 Preliminaries 

It is assumed that the readers are familiar with the basics of Membrane Comput- 
ing^] and formal language theory. We only briefly mention notions and notations 
which will be useful throughout the paper. Let V be an alphabet, V* is the free 
monoid over V with respect to concatenation and the identity element A (the 
empty string). The set of all non-empty strings over V is denoted as V + so 
V + — V* — {A}. We call V a singleton if V = {a} and simply write a* and 
a + instead of {a*} and {a + }. The length of a string w £ V* is denoted by 
\w\. If a is a symbol in V, a Q — A. A language L C V* is regular if there is 
a regular expression E over V such that L(E) = L. A regular expression over 
an alphabet V is constructed starting from A and the symbols of V using the 
operations union, concatenation, and +, using parentheses when necessary to 
specify the order of operations. Specifically, (i) A and each a £ V are regular 
expressions, (ii) if Ei and E 2 are regular expressions over V then (Ei U E%), 
E\E%, and E^ are regular expressions over V, and (tit) nothing else is a reg- 
ular expression over V. With each expression E we associate a language L(E) 
defined in the following way: (i) L(X) = {A} and L(a) = {a} for all a £ V, (ii) 

2 A massively parallel SNP system simulator in [213] is based on the matrix represen- 
tation in order to perform simulations. 

3 a good introduction is [10| with recent results and information in the P systems 
webpage at http : / /ppage . psystems . eu/ and a recent handbook in [12] 



L(E 1 U E 2 ) = L(E X ) U L(E 2 ), L(E 1 E 2 ) = L(E 1 )L(E 2 ), and L(E+) = L(E t )+, 
for all regular expressions E 1: E 2 over V . Unnecessary parentheses are omitted 
when writing regular expressions, and E + U {A} is written as E* . 
We define an SNP system of a finite degree m > 1 as follows: 

II = (0,ai,...,a m ,syn), 

where: 

1.0 = {a} is the singleton alphabet (a is called spike). 

2. ai, . . . , a m are neurons of the form <x; = [on, Ri),l < i < m, where: 

(a) ai > is an integer representing the initial number of spikes in <7j 

(b) Ri is a finite set of rules of the general form 

E/a c -» a b ;d 

where E is a regular expression over O, c > 1, 6 > 0, with c > b. 

3. syn C {1, 2, . . . , m} x {1,2,..., to}, (i, i) ^ syn for 1 < i < m, are synapses 
between neurons. 

A spiking rule is where b > 1. A forgetting rule is a rule where 6 = is 
written as i?/et c — ► A. If L(-E) = {a c } then spiking and forgetting rules are 
simply written as a c —> a b and a c —> A, respectively. Applications of rules are 
as follows: if neuron ai contains k spikes, a k £ L(E) and k > c, then the rule 
E/a c — ¥ a b € Ri is enabled and the rule can be fired or applied. If b > 1, the 
application of this rule removes c spikes from ai , so that only k — c spikes remain 
in ai. The neuron fires b number of spikes to every Gj such that (i,j) S syn. If 
the delay d = 0, the b number of spikes are sent immediately i.e. in the same 
time step as the application of the rule. If d > 1 and the rule with delay was 
applied at time t, then the spikes are sent at time t + d. From time t to t + d — 1 
the neuron is said to be dosecQ and cannot receive spikes. Any spikes sent to the 
neuron when the neuron is closed are lost or removed from the system. At time 
t+d the neuron becomes open and can then receive spikes again. The neuron can 
then apply another rule at time t + d+1. If 6 = then no spikes are produced. 
SNP systems assume a global clock, so the application of rules and the sending 
of spikes by neurons are all synchronized. 

A configuration of the system can be denoted as (rit/ti, ■ ■ ., n m /t m ), where 
each element of the vector is the configuration of a neuron Oi , with spikes and 
is open after t{ > steps. At time step 0, since no rules, with or without delay, 
are yet to be applied, tj = 0,Vi € {1, ...,to}. U increases when ai applies a rule 
with delay. 

To illustrate, let us have an snp system with three neurons, a\, a 2l and 03 
with synapses (1, 2), (2, 3). Only o\ has some number of spikes, say x. o\ has a 
rule a + / a — > a; 2, a 2 has a rule a — > a, and 173 is a sink neuron. At a time step k 
where no rules are yet to be applied, the configuration of this system would be 
(x/0, 0/0, 0/0). At time step k + 1, we apply the rule in a\\ the configuration 



4 this corresponds to the refractory period of the neuron in biology 



Fig. 1. An example of an SNP system where the source neuron is <r\ and sink 
neuron is 05. Neuron 174 will only spike once it accumulates two spikes, one each 
from (72 and 173. 

would then become (x- 1/2, 0/0, 0/0). At time step k + 2, (x-1/1, 0/0, 0/0). 
At time step k + 3, o\ is open again and will now fire one spike to 02, the 
configuration would now be (x — 1/0, 1/0, 0/0). At time step k + 4, o\ can 
again use its rule, and now 02 can also apply its rule; the configuration would 
then become (x - 2/2, 0/0, 1/0). 

The initial number of spikes of o~i is denoted as on, on > 0. The initial configu- 
ration therefore is (a±/0, . . . , a m /0). A computation is a sequence of transitions 
from an initial configuration. A computation may halt (no more rules can be 
applied for a given configuration) or not. 

For SNP systems in this worl^Jwe have only one source neuron (a neuron 
without any incoming synapses to it) and one sink neuron (a neuron without any 
outgoing synapses from it). The initial configuration will always be a single spike 
found only at the source neuron. The objective is to route this spike from the 
source to the sink neuron. The routes will involve paths in the system, where a 
path consists of at least two neurons <7j, o~j such that <E syn. Using this idea 
of a path, we can have four basic routing constructs: (1) sequential, where <Ji+i 
spikes after o~i and (i,i + l) G syn, (2) iteration, where at least two neurons spike 
multiple (possibly an infinite) number of times and a loop is formed e.g. synapses 
and exist between neurons <7j and Uj, (3) join, where spikes from at 
least two input neurons o~ m , a n are sent to a neuron o~i, where (m, i), (n, i) £ syn, 
so that Oi produces a spike to at least one output neuron aj, and € syn, 
(4) split, where a spike from Uj is sent to at least two output neurons <7fc and 07 
and (j, k), (j, I) € syn. An example of an SNP system that we deal with in this 
work is shown in Fig. [T] 

SNP systems considered in this work are those with neurons having exactly 
one rule only. For SNP systems with delays, notice that if there exists a sequential 
path from <7; (with delay dl) to o~j (with delay d2) if dl < d2 and o~i spikes more 
than once, it is possible for the spikes from o~i to be lost. This is due to the 
possibility that o~j may still be closed when spikes from o~i arrive. We therefore 
avoid lost spikes by assuming dl > d2 whenever a sequential path from o~i to <7j 

5 more information on SNP systems having input and output neurons as well as non- 
deterministic SNP systems (we focus on deterministic systems in this work) can be 
found in [11] and [12] for example. 



exists, and leave the cases where dl < d2 as an open problem. Given an SNP 
system with delays 77 and an SNP system without delay 77, by simulation in 
this work we mean the following: (a) the arrival time t of spikes arrive at sink 
neurons of 77 is the the same time t for the arrival of spikes at the sink neurons 
of 77 or offset by some constant integer k, (b) the number of spikes that arrive 
in the sink neurons of 77 is the same for the sink neurons in 77 or a factor of the 
delay of 77. We define the total runtime of 77 and 77 as the total time needed 
for a spike to arrive to the sink neurons from the source neurons. 



3 Main Results 

Now we present our results in simulating SNP systems with delays using SNP 
systems without delays. Once again we denote SNP systems with delays as 77 
and those without delays as 77, subject to the restrictions and assumptions 
mentioned in the previous section. 

Definition 1. Given an SNP system, II = (0,ax, . . . ,o~ m , syn), a routing of 
77, 77' is defined as 

77' = (0,E',syn'), 

where: 

1. syn' C syn 

2. E is a subset of the set of neurons of II. E' = {0, <r|(er, a) G syn'}. 

Lemma 1. Given an SNP with delay, II that contains a sequential routing, 
there exists an SNP without delay, U that contains a routing that can simulate 
the sequential routing of II . 

Proof. We refer to Fig. [2] for illustrations and designate empty neurons as sink 
neurons. Let 77 be the system in Fig. [2]Ja). Total runtime for 77 is d steps, where 
final spike count for er 12 is ai2 = 1- Now let 77 be Fig. [2^b) that simulates 77. 
The initial spike count for 021 is 0:21 = 1 + d. Total runtime for 77 is 1 + d steps, 
and 022 receives 1 + d number of spikes. Final spike count for 023 is 023 = 1. 

Now let 77 be the system in Fig.J^c) and 77 be the system in Fig.^d). 77 now 
has two delays, dl and d2, and the total runtime for 77 is 1 + g?1 + g?2. For 77 (with 
initial spikes a = (1 + dl)(l + d2j), the total runtime is (1 + dl)(l + d2) +3 - 1. 
The always open neurons of 77 except for 044 each introduce one time delay i.e. 
3, minus one time delay since 044 does not spike. The final spike count for 033 
is a = 1, which is also the final spike count for 044. It can be easily shown that 
addition of neurons without delayed rules in either 77 or 77 contribute to a delay 
of one time step (per neuron) being added to the total runtimes for both systems 
(see Fig. [3). □ 

The total runtime for 77 is one additional time step from the total runtime 
of 77. The additional time step for 77 comes from delay due to the first spiking 
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Fig. 2. Sequential routing, with single and multiple delays. 
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Fig. 3. Sequential routing with additional neurons without delays. 



Lemma 2. Given an SNP with delay, U that contains an iterative routing, 
there exists an snp without delay, U that contains a routing that can simulate 
the iterative routing of II . 



Proof. We refer to Fig. [4] and [5] for illustrations. Let 77 be the system in Fig. 
Qa) and II be the system in Fig. |4|b). Notice that o\\ and u\2 continuously 
replenish the spike of one another, forming an infinite loop so that a±3 will keep 
on accumulating spikes. The total runtime is 1 + d steps. Neurons a 2 i and CT22 
similarly replenish the spike of one another in an infinite loop. The total runtime 
for II is 2 + d, with one additional delay due to CT23. Due to the operation of 023, 
both <7i3 and 024 obtain one final spike each. It can be easily shown that if the 
delay of 77 is at an instead of 012, the 77 in Fig. |4pt>) can still simulate 77. 

For iterations with more than one delay, let 77 be the system in Fig. [SJa) 
and 77 be the system in Fig. [f^b). The total runtime of 77 is dl + d2 and 
(1 + dl) (1 + (72) + 3 — 1 for 77. The final spike count for 77 is 0113 = 1 which is 
the same as the final spike count for 77 (due to the operation of 0-24). □ 
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Fig. 4. Iterative routing. Note that the 77 in (b) that simulates 77 in (a) is still 
the same even if an has the delay instead of a 12 



Notice that the total runtime for 77 for iterative routing with multiple delays 
is the same as the total runtime for 77 for sequential routing with multiple delays. 
Again, every additional neuron without delay adds one time delay to the total 
runtime of 77 and 77. 

Lemma 3. Given an SNP with delay, II that contains a split routing, there 
exists an SNP without delay, II that contains a routing that can simulate the 
split routing of II . 

Proof. We refer to Fig. [6] and [7] for illustrations. Let 77 be the system in Fig. 
[HJa) and 77 be the system in Fig. [6|b) and we first consider a split where the 
neuron on that performs the split is the one that has a delay. The total runtime 
for 77 is 1 + d, with a final spike count of 0:13 = a\4 = 1. The total runtime for 
77 is also 1 + d with a similar final spike count of 0:23 = «24 = L 

Now let 77 be the system in Fig. |7|a) and 77 be the system in Fig. [7^b) and 
we consider a split where the delay is in one of the receiving child neurons (013). 
In this case the sink neurons do not receive a spike at the same time. For 77, if 
till spikes at time t, 014 receives a spike at time t + 1 while 015 receives a spike at 
time t+l + d. The total runtime therefore (time when both sink neurons obtain 
one spike) is 1 + d. For 77 if 021 first spikes at time t then 024 first receives a 
spike at t+1. The total runtime is 1 + d where a 2 5 = 1 and a 2 4 = 1 + d. □ 

Note that since 024 accumulates 1 + d spikes (from 021), if the split SNP 
system 77 is a subsystem of a larger SNP system, an additional neuron erj with 
a rule a 1+d — > a needs to be added in order to return the spike number back to 
oncEl 

Lemma 4. Given an SNP with delay, II that contains a join routing, there 
exists an SNP without delay, II that contains a routing that can simulate the 
join routing of II. 

Proof. We refer to Fig. [8] and [9]f or illustrations. Let 77 be the system in Fig.|8ja) 
and 77 be the system in Fig.]8[b) so that we first we consider a join where the 



this additional neuron once again adds one time step of delay. 
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Fig. 5. Iteration with more than one delay. 



delay is in a neuron <Ti2 before the neuron an that performs the join. The total 
runtime for 77 is 1 + d since CT13 will wait for the spike from parent neuron <Ji2 
to arrive before sending a spike to CT14. The total runtime for 77 is 2 + d because 
of the additional source neuron <7 2 ij since 1724 will wait for the spike from er 2 3. 
For II and U the final spike count is 1. 

Now let II be the system in Fig. [9^a.) and 77 be the system in Fig. j9^b) and 
we consider next a join where the delay is in the neuron (T13 that performs the 
join. Total runtime for 77 is 1 + d which is also the total runtime for 77. The 
final spike count for both systems is 1. □ 
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Fig. 6. A split where the neuron performing the split has a delay. 
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Fig. 7. A split where a child neuron has a delay. 
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Fig. 8. A join where a parent neuron has a delay. 



Theorem 1. Let LI be an SNP system with delays that follows the restrictions 
and assumptions given in Section^ and contains routings limited to the following: 

1. sequential 

2. iterative 

3. split 
4- join 

Then, there exists an SNP system without delay, LL , that simulates LL . 
Proof. Proof follows from Lemma 1, 2, 3, and 4. □ 



4 Final Remarks 

Since LJ has to simulate 77, it is expected that either additional neurons, regular 
expressions, or spikes are to be added to LI in order to simulate the delay(s) 



a 

a + /a — » a 



i-V 



(a) 



12, 



(a 2 ) + /a — > 



14 









2 V 



(b) 



a + /a — > a 



,23 



22 ^ a l+d2 



a + /a — » a 



(a 



24 



Fig. 9. A join where the neuron performing the join has a delay. 



and final configuration of 77 (with some offset or spike count multiple of II). We 
have shown in this work how to the routing constructs of II, although we could 
perhaps do better i.e. less neurons, less initial spikes, and so on. 

Conjecture 1. For the case a = a k where k > 1 is an integer, SNP systems 
without delays can still simulate SNP systems with delays. 

Conjecture 2. It is possible to make the time of spiking of SNP systems with 
delays exactly coincide with the time of spiking of SNP systems without delays 
(or at least to lessen the neurons or time difference as presented in this work). 

We did not consider a split where the delays of the child neurons are not 
equal, since Lemma [3] only considers a split where only one of the child neurons 
have a delay (also, it is easy to show that if both child neurons of a split have the 
same delay d then the H in Lemma [3] still holds). As mentioned earlier, we also 
leave as an open problem on how to simulate SNP systems with delays where 
dl < d2 for sequential routing. 
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